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ABSTRACT 


Rayleigh  wave  data  were  recorded  by  a  single  vertical 

high-gain  seismograph  located  in  north  Texas.  It  was  found 

that  Rayleigh  waves  propagating  in  a  continental  patn  over 

the  pole  to  this  station  exhibited  very  stable  dispersion 

characteristics,  particularly  in  the  period  range  35-75  secs: 

at  shorter  periods  an  average  dispersion  curve  could  be 

astimated.  A  composite  dispersion  curve  was  used  to  develop 

"chirp"  (or  matched)  filters,  the  form  of  which  depended  on 

the  distance  between  the  seismograph  and  the  event.  Sucn 

filters  are  shown  to  be  very  efficient  in  improving  the 

signal-to-noise  ratio  of  Rayleigh  waves  emanating  from  events 

in  the  Sino-Soviet  region  and  can  be  used  as  a  means  of 

separating  "mixed"  Rayleigh  waves.  For  the  path  over  the 

pole  (or  "polar  waveguide")  the  detection  threshold  is  estimated 

to  be  at  M  =3.5. 
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Rayieigh  wave  da. a  were  recorded  by  a  single  vertical  high  gain  seismograph  located 
in  north  Texas.  It  was  found  that  Rayleigh  waves  propagating  in  a  continental  path 
over  the  pole  to  this  station  exhibited  very  stable  dispersion  characteristics, 
particularly  in  the  period  range  35-75  secs:  at  shorter  periods  an  average 
ispersion  curve  could  be  estimated.  A  composite  dispersion  curve  was  used  to 
develop  chirp  (or  matched)  filters,  the  form  of  which  depended  on  the  distance 
between  the  seismograph  and  the  event.  Such  filters  are  shown  to  be  very  efficient 
xn  improving  the  signal- to-noise  ratio  of  Rayleigh  waves  emanating  from  events  in  the 
ino-boviet  region  and  can  be  used  as  a  means  of  separating  "mixed"  Rayleigh  waves. 

F jr  the  path  over  the.  pole  (or  "polar  waveguide")  the  detection  threshold  is 
estimated  to  be  at  M  =3.5. 
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I  12.  SPONSORING  MILITARY  ACTIVITY 


INTRODUCTION 

Ewing  and  Press  (1952)  described  the  classical  "peak- 
and-trough"  method  of  determining  a  fundamental  dispersion 
curve  from  a  single  surface  wave  coda.  Sato  (1955)  applied 
Fourier  analysis  to  the  calculation  of  dispersion  curves. 

These  techniques  have  been  improved  upon  (Block  and  Hales, 
(1968),  Landisman,  et  al  (1969)  and  Dziewonski,  et  al  (1969) ); 
and  an  overall  processing  system  was  devised  combining  Fourier 
analysis  and  time-varying  filters,  enabling  the  calculation  of, 
not  only  the  fundamental  phase  and  group  velocity  curves,  but 
also  the  curves  for  the  higher  modes.  In  particular,  once  a 
group  velocity  curve  has  been  obtained  it  is  possible  to 
design  a  time  varying  filter  that  extracts  a  single  dispersed 
mode  and  effects  a  significant  improvement  in  the  signal-to- 
noise  ratio.  The  filtered  output  is  then  an  estimate  of  the 
true  signal.  Dziewonski,  et  al  (1968)  summarizes  these 
methods  and  it  is  evident  that  they  extract  the  maximum  amount 
of  information  from  the  surface  wave  signal.  Nevertheless, 
the  methods  are  both  time  consuming  and  computationally 
difficult.  But  if  we  restrict  the  investigation  to  the 
fundamental  mode  of  surface  wave  propagation,  it  becomes 


possible  to  develop  a  fast,  simple  digital  technique  for 
the  analysis. 

In  this  paper  we  describe  a  scheme  which  estimates  the 
fundamental  group  velocity  dispersion  curve  from  a  seismo¬ 
gram,  and  then  constructs  a  "matched"  filter  based  on  the 
estimated- curve.  This  matched,  or  chirp,  filter  can  then 
be  used  to  operate  on  signals  with  low  signal-to-noise  ratios 


THE  SEISMOGRAMS 


During  the  period  March  through  August  1970  a  three- 
component  seismograph  system  was  operated  at  a  depth  of 
180  m  in  the  Morton  Salt  Company  mine  at  Grand  Saline,  Texas 
(Sorrells  et  al. ,  1971).  The  seismometers  were  sealed  in 
such  a  manner  that  the  masses  were  not  subject  to  buoyancy 
effects  from  atmospheric  pressure  changes.  The  systems  were 
of  the  "advanced"  long  period  type  which  are  operated  with 
a  peak  magnification  at  a  period  of  50  secs.  Sorrells  et 
al .  (1971)  have  shown  that  if  such  a  system  is  operated  at 
a  site  which  is  remote  from  seismic  noise  generated  by  at¬ 
mospheric  processes,  the  spectra  of  the  long  period  noise 
exhibits  a  minimum  in  the  period  range  20-60  sec.  Careful 
installation  of  these  seismographs  at  a  depth  of  about  180  m 
enabled  the  vertical  seismograph,  the  only  one  used  in  the 
present  study,  to  be  operated  at  a  magnification  of  200,000. 
The  data  were  recorded  on  a  digital  acquisition  system  de¬ 
scribed  by  Herrin  and  McDonald  (1971)  after  sampling  at  a 
rate  of  1  per  second. 

The  recorded  seismograms  were  associated  with  source 
parameters  and  body  wave  magnitudes  provided  in  Preliminary 
Determination  of  Epicenters  ( PDE)  issued  by  NOAA. 
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DISPERSION  CURVE  ESTIMATION 


A  computer  program  was  written  which  enabled  analog 
representations  of  digital  tine  series  to  be  displayed  on 
an  oscilloscope.  Selected  seismograms  were  then  plotted, 
using  a  Cal-comp  plotter.  It  was  decided  to  limit  the  ob¬ 
servations  to  Rayleigh  waves  recorded  at  teleseismic  distance 
which  had  traversed  a  predominantly  continental  path;  thus 
restricting  'the  range  of  azimuths  to  about  15°  east  and  west 
of  the  pole,  and  the  range  of  epicentral  distances  to  55°-130°. 

Using  the  plot  of  the  seimogram  an  estimation  was  made 
of  the  start  time  (Ts)  and  the  length,  in  seconds,  (window) 
of  the  Rayleigh  signal.  A  computer  program  was  written  which 

searched  this  window,  starting  at  T  ,  for  zero  crossings,  in 

s 

a  manner  similar  to  the  "peak  and  trough"  tech,  ique  of  Ewing 
and  Press,  (1952).  The  time  of  the  "actual"  zero  crossing 
was  calculated  by  triangulation  between  the  last  point  immedi¬ 
ately  prior  to  a  zero  crossing,  and  the  first  point  immediately 
after  it.  The  initial  zero  crossing  time,  (T  )  was  stored,  and 
at  the  next  zero  crossing,  (T^  the  first  period,  (Px)  was 

computed.  T  and  P,  were  then  stored  and  the  seach  continued 
1  x  A 

until  the  window  was  exhausted.  For  i  =  1,  2,  ...,  N,  the 
group  velocites  and  the  periods  were  computed  from 
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I 


VU=A  ^  PL=  2  (TrT.,)^ 

T_ri. 

where  N  =  the  number  of  zero  crossings  in  the  window, 

A  =  the  epicentral  distance  from  event  to  station  in 
km,  and 

TT^  =  the  observed  travel  time  from  source  to  the  ith 
zero  crossing. 

Since  the  actual  zero  crossing  occurred  within  a  one  second 

interval  any  variation  must  be  less  than  one  second.  As  a 

conservative  approximation,  we  assumed  the  error  to  be  uniform; 

2 

thus  the  variance  was  1/12  sec  . 

Now  TT.  =  T.  -  0 
l  i 

where  0  is  the  origin  time  of  the  event.  From  Herrin  et  al. 
(1968)  a  reasonable  origin  time  variance  is  0.81  to  1.00  sec2. 
Since  and  0  are  independent  observations  their  variances 
may  be  added.  We  thus  take  1.00  sec2  as  a  reasonable  variance 


for  each  T1\  .  From  Tucker  et  al.  (1968)  and  Herrin  et  al.  (1968) 
a  reasonable  variance  in  A  is  900  km2.  We  then  have 


»i 

TTL+ 


mMtt 


where 


^S.  =  the  true  distance  (in  km) 


TT.  =  the  true  travel  time. 
1  ' 


el  =  the  error  in  epicentral  distance  (  0"1=3O  km), 


e2  =  the  error  in  travel  time  (  C"2=  1  sec) 


Assuming  an  event  with  ^  =  6000  km  and  a  maximum  group 


velocity  of  4.00  km/sec  (Oliver,  1962)  then 


V  _  (6000+  e  )  / 

1  /CIS00+  ea)  ’ 


Since  the  variance  of  e2  is  small  relative  to  1500  we  can 


ignore  e^  and  write 


Vt=  4.o+  e^/isoo 
”  4".  o  +  £  ^  , 


where 


CT  =  0  0004-  sec1 . 

S1 


For  a  lower  period  cutoff  of  20  sec  (or  a  velocity  of  3.00 


km/sec)  we  get 


where 


Vn*  (GOOO+  ej^000 
=  3  o  +  e„  / 200  o 
=  3-o  + 


CT  a  0-000  22  S  set1. 


Thus  as  a  conservative 


where 


1 


approximation  we  can  assume 
=  True  velocity  i  +  £ 

=  0.0004  sec2. 


Let  us  now  consider  the  variation  in  the 


We  have 


P. 

1 


wnere  o 

»  O  OS33 
L 

Since  the  minimum  period  is  assumed  to  be  20  sec  the  minimum 
time  difference  is  of  the  order  of  10  sec  and  we  may  assume 
that  the  errors  in  T^  and  T^_^  are  independent.  Then  the 
variance  of  a  Pi  is 

Cp=4Co.°?33f  °  0S33)  =  0.GGS4  iec1. 

'i 

Thus  the  variance  in  a  velocity  observation  is  three  orders 
of  magnitude  less  than  that  of  a  period  observation.  We 
therefore  ignored  the  error  in  the  velocity  computations  and 
assumed  that  the  velocities  were  known  without  error.  Now 
assuming  that  all  the  error  lies  in  the  period  computations, 
it  is  possible  to  consider  a  group  velocity  curve  as  a 
regression  of  period  on  group  velocity. 

Figure  2  shows  the  data  points  obtained  from  dispersed 
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Rayleigh  waves  for  three  events  located  in  the  Arctic  (Table  1) 
at  similar  epicentral  distances  and  similar  azimuths.  Data 
from  an  event  in  Sinkiang  Province  along  a  similar  azimuth, 
but  at  a  greater  distance,  (Table  1)  are  shown  in  Figure  3. 

A  comparison  of  the  data  in  Figures  2  and  3  in  the  velocity 
range  3.50  to  3.90  km/sec  (periods  approximately  35  and  75  sec 
respectively)  indicates  that  there  is  no  statistical  difference 
between  the  true  dispersion  curve  for  each  event.  The  observed 
scatter  could  result  from  statistical  fluctuations  in  the 
estimates  or  real  differences  in  the  velocity-period  relation¬ 
ship  along  the  propagation  path.  Thus  it  appears  that  in  the 
velocity  range  3.50  to  3.90  km/sec  there  is  a  unique  disper¬ 
sion  curve  independent  of  path. 

However,  in  the  velocity  range  2.90  -  3.5  km/sec  (periods 
of  20-35  sec) ,  there  are  indications  that  real  differences  may 
exist  between  the  true  dispersion  curves  in  Figures  2  and  3. 

In  this  velocity  range  we  estimated,  for  the  purposes  of  this 
study,  an  average  dispersion  curve. 

Using  group  velocity  as  the  independent  variable  the 
data  for  each  event  were  grouped  into  cells  of  width  0.5  km/sec. 
Within  each  cell  the  data  were  ordered  according  to  period  value 
end  the  sample  median  was  computed.  The  cell  midpoint 


velocities  and  median  periods  were  stored  as  data  pairs. 

The  sample  median  was  chosen  as  a  measure  of  central  tendency 
because  it  is  more  efficient  than  the  sample  mean  in  the 
presence  of  outliers  (Dixon,  1953) 

In  the  velocity  range  3.50  to  4.00  km/sec  the  sample 
median  periods  are  observations  on  a  single  curve;  therefore, 
adjacent  points  are  related  and  a  smoothing  operator  should 
improve  the  results.  Furthermore,  in  the  velocity  range 
2.90  to  3.50  km/sec  only  a  representative  smooth  curve  was 
desired,  and  a  three  point  moving  average  operator  was 
applied  to  the  sample  period  values.  The  operator  is  given  by 

Period  (1)  =  x (1) 

Period  (Num)  =  1/2  •  X(Num-l)  +  1/2  .  X(Num) 


Period  (i)  =  1/4  •  X(i-i)  +  1/2  .  x(i) 

+1/4  *  X  (i+1) 


where 


1-2,  . . . ,  Num-1 

1  =  the  index  of  the  cell  with  the  largest 
velocity 

Num  =  the  index  of  the  cell  with  the  smallest 
velocity 

X(i)  =  the  sample  median  period  of  cell  i 
Period  (i)  =  the  smoothed  period  of  cell  i. 
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The  resulting  dispersion  curves  for  the  four  events 
(Figures  2  and  3)  are  shown  in  Figure  4.  The  smoothed  curves 
show  the  same  results  as  the  original  zero  crossing  data. 

That  is,  the  curves  are  not  statistically  different  in  the 
velocity  range  3.50  to  3.90  km/sec,  and,  in  the  range  2.90 
to  3.50  km/sec,  the  smoothed  curves  indicate  that  the  actual 
paths  do  not  produce  markedly  different  dispersions  curves 
at  these  short  periods. 

In  view  of  the  similarity  of  the  curves  we  combined 
the  data  from  the  four  events  and  estimated  a  composite 
group  velocity  dispersion  curve  shown  in  Figure  5.  Also  shown 
in  this  figure  is  the  continental  dispersion  curve  compiled 
by  Oliver  (1961);  the  long  dashed  lines  indicate  bounds  of 
the  expected  scatter  under  normal  observational  conditions. 

The  composite  curve  can  be  seen  to  fall  within  the  error 
bounds.  Furthermore,  the  composite  curve  is  in  good  agree¬ 
ment  with  the  results  of  Ewing  and  Press  (1956)  who  reported 
a  group  velocity  of  about  3.9  km/sec  at  a  period  of  75  sec. 

The  smoothed  data  points  are  given  in  tabular  form  in  Table  2. 
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CHIRP  FILTER  CONSTRUCTION 

We  have  shown  that,  for  our  seismograph  station  in 
north  Texas,  we  can  assume  a  single  group  velocity  dispersion 
curve  for  the  continental  path  over  the  pole.  We  can  then 
design  "matched",  or  chirp,  filters  to  aid  in  the  detection 
of  surface  waves  from  small  earthquakes  traversing  the  same 
path. 

A  program  was  written  which  produced  a  sine  wave  that 
decreased  in  period  with  time,  the  dispersion  being  a  function 
of  epicentral  distance.  Such  waveforms  have  been  given  the 
name  "chirp"  in  radar  technology.  In  this  case  the  informa¬ 
tion  used  to  distance  -  function  were  from  Table  2.  The 
following  assumptions  were  made: 

(i)  the  first  arriving  energy  is  at  the  maximum  period 
observed,  anc 

(ii)  between  the  data  points  of  Table  2  the  dispersion 
curve  is  linear. 

With  these  assumptions  we  obtained,  by  simple  algebraic 
manipulations, 

TV  %  W  To=TTl“  V* 

VL 


and  for  the  ith  linear  segment 

-W.Jfc)  -  2tt  CtfT^C1/^.,-  Vtt  J _ 

[P.?  V  (p.-,  /IT,-  n.  /TT^Mt^T.)] 

“here  CTL.rT.U  ±  , 

=  frequency  in  radians  for  the  t  th  time, 

=  the  i  th  period  in  sec, 

TT^ >  =  the  i  th  travel  time  in  sec, 

i  =  1 ,  2 ,  . . . ,  Num . 

The  dispersed  sine  wave  computations  were  given  by 
AMP  (k)  =  sin  k  +  f  i']  , 

where  k  =  1,  2,  . . . ,  Ni 

=  the  number  of  data  points  in  the  i  th  segment 
Pi  =  phase  shift  of  the  i  th  segment 
AMP (k)  =  the  dispersed  sine  wave  amplitude. 

For  k  =  1  the  segment  was  extrapolated  to  include  time  points 
back  to  Tq.  For  all  other  values  of  k  only  time  points  be¬ 
tween  the  TT^  were  included  in  a  given  segment. 

The  p  were  as  follows 

Pl=° 

Pirn  (  +  Pu-|  » 
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for  1  =  2,  3,  . ..,  Num,  whtre  t  is  the  largest  t  value  in 

max 

the  i  th  segment.  Also 

fc=  k+  TT 


so  that 


•fc 

W16*. 


In  each  linear  segment  the  frequency  was  a  function  of 
time  given  by  TATV  (t).  The  origin  of  the  sine  wave  was  taken 
as  zero,  which  gave  the  time  shift  between  t  and  k.  The 
phase  shift  p  ^  ensured  a  smooth  transition  from  segment  to 
segment . 

In  Figure  6  we  show  the  dispersed  sine  wave  for  one 
of  the  events  used  in  calculating  the  composite  dispersion 
curve  (Figure  5) .  These  sinusoids  indicate  the  effects  of 
differences  in  epicentral  distance;  for  the  Sinkiang  Province 
event  (A  =  107.6°)  the  first  zero  crossing  occurs  at  about 
72  sec  and  for  the  other  three  (  &  =  approximately  60°)  the 
first  zero  crossing  occurs  at  about  68  sec.  Similar  period 
differences  are  evident  throughout  the  curves. 

Convolution  of  these  dispersed  sine  waves,  or  chirp 
filters,  with  their  respective  signals  should  then  produce 
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psuedo -autocorrelation  functions.  The  results  are  not  true 

\ 

autocorrelations  because  the  sine  waves  have  not  been 
amplitude  modulated;  nor  is  the  "true"  dispersion  given 
exactly  by  the  composite  dispersion  curve. 


APPLICATION  OF  THE  CHIRP  FILTER 

\ 

The  chirp  filter  constructed  from  the  dispersion  curve 
in  Figure  5  was  applied  to  four  different  sets  of  seismograms. 
First,  it  was  convolved  with  the  original  seismograms  used  to 
construct  the  dispersion  curve;  second,  with  other  events 
remote  from  the  first  four;  third,  as  detector  of  surface 
waves  with  a  low  signal-to-noise  ratio;  and  fourth,  as  a 
means  of  separating  "mixed"  Rayleigh  waves. 

In  order  to  remove  very  long  period  energy  each  seismo¬ 
gram  was  filtered  prior  to  being  convolved  with  the  chirp 
filter.  The  high-pass  digital  filter  used  had  a  corner  at 
50  sec  and  fell  at  36  db/octave  to  a  minimum  at  a  period  of 
100  sec. 

Figure  7  shows  the  application  of  the  chirp  filter  to 
the  time  series  for  each  of  the  events  used  in  calculating  the 
composite  dispersion  curve.  In  each  example  the  lower  trace 
shows  the  filter  output,  correctly  aligned  with  the  input 
time  series.  The  auto-correlation-like  nature  of  the  output 
should  be  noted;  the  maximum  peak  of  which  should  occur  at  the 
beginning  of  the  Rayleigh  wave  energy.  In  each  of  the  figures 
the  data  were  plotted  so  that  the  largest  excursion  was  full- 
scale;  therefore,  the  improvement  in  signal-to-noise  ratio 


can  be  estimated  by  examining  tne  reduction  in  noise. 

Three  other  events  were  selected  (table  3)  with  similar 
continental  polar  surface  wave  paths  to  Texas;  however, 
their  epicenters  were  remote  from  those  given  in  table  1. 

The  results  of  convolving  these  seismograms  with  the  same 
chirp  filter  used  previously  are  shown  in  figure  8.  The 
outputs  clearly  show  that  the  filter  is  effective  for  widely 
separated  sources. 

Figure  9  shows  the  results  of  convolving  the  chirp 
filter  with  seismograms  in  which  the  surface  waves  are 
barely  discernible.  This  figure  shows  a  seismogram  from  a 

small  event  (  rw  *  4.7)  at  a  depth  of  110  km  in  the  Rindu- 

Ir 

Kush;  convolution  clearly  emphasized  the  start  of  the  surface 
waves. 

Mixed  Rayleigh  waves  resulting  either  from  multiple 
sources  close  together  or  from  reflections,  or  from  multi- 
pathing  are  very  difficult  to  separate,  even  when  data  from 
large  arrays  are  avialable.  However,  if  the  arrivals  are 
separated  by  times  greater  than  about  50  sec,  the  chirp 
filter  technique  can  be  useful  in  separating  the  arrival  . 
Figure  10(a)  (lower  trace)  shows  the  filtered  output  for  a 
small  event  from  the  region  of  the  USSR-Mongolia  border; 
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two  and,  possibly,  three  Rayleigh  wave  arrivals  can  be  seen 
on  the  filtered  trace.  The  first  and  largest  of  these 
arrivals  occurs  precisely  at  the  expected  arrival  time  of 
the  fundamental  Rayleigh  wave  for  this  event.  A  fairly 
large  earthquake  (tn  =  5.9)  sparked  an  earthquake  swarm  from 
this  region  of  the  USSR-Mongolia  border  in  May  1970. 

Another  earthquake  from  the  swarm  is  shown  in  figure  10(b); 
again  convolution  detected  two  arrivals.  All  available 
seismograms  of  earthquakes  in  the  swarm  were  processed  but 
only  three  showed  multiple  detections  of  surface  waves. 


DISCUSSION 


We  have  described  above  a  new  approach  to  the  "peak 
and  trough"  method  of  surface  wave  analysis  which  has  been 
used  to  determine  group  velocity  dispersion  curves.  The 
machine  measurements  proved  to  be  at  least  as  good  as  those 
of  an  experienced  analyst,  and  could  be  produced  much  faster. 

A  further  advantage  was  gained  in  the  present  analysis 
from  the  special  properties  of  the  seismograph  system. 

Many  previous  surface  wave  studies  have  used  data  from  the 
World-Wide  Standard  Seismograph  Network  recordings  in  which 
the  peak  period  of  the  calibration  response  is  at  20  sec. 

In  our  systems,  the  peak  response  proved  to  lie  within  the 
"stable "portion  of  the  group  velocity  dispersion  curve. 

This  type  of  analysis,  while  being  quick  and  computa¬ 
tionally  easy,  nevertheless  has  its  problems.  Any  non-least- 
time  arrivals  will  impair  the  estimation  of  the  dispersion 
curve;  the  technique  does  not  give  as  much  information  as 
that  of  Dziewonski  et  al .  (1969),  which  places  a  restriction 
on  the  use  of  the  method. 

Nevertheless,  by  carefully  selecting  events  with  a 
large  signal-to-noise  ratio,  and  with  an  absence  of  either 
higher  modes  or  very  long  period  surface  waves,  it  has  been 
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shown  that  a  single  fundamental  group  velocity  dispersion 
curve  can  be  calculated  for  a  continental  path  over  the 
pole  to  the  central  United  States.  This  claim  can  be  further 
substantiated  by  the  fact  that  a  chirp  filter,  constructed 
from  this  fundamental  dispersion  curve,  proved  effective 
for  events  from  a  region  extending  from  southern  USSR  and 
northern  China  to  the  Arctic  Ocean. 

In  addition  to  the  events  mentioned  in  the  previous 
sections  of  this  paper,  seismograms  from  the  other  events 
listed  in  table  3  were  processed.  In  some  cases  the  auto¬ 
correlations  displayed  low  signal-to-noise  ratios,  but  in 
every  search  for  an  event  reported  by  NOAA,  an  event  was 
detected?  unreported  events  were  also  detected.  Thus  further 
study  is  required  to  determine  operational  detection  thres¬ 
holds  and  false  alarm  rates. 


CONCLUSIONS 


In  view  of  the  results  reported  here  it  might  be  appro¬ 
priate  to  suggest  a  reason  for  the  stability  of  the  group 
velocity  dispersion  curve,  particularly  at  the  longer  periods, 
on  a  continental  path  over  the  pole. 

It  has  been  shown  in  recent  years  (e.g.  Roy  et  al.  1972) 
that  the  western  United  States,  known  structurally  as  the  Basin 
and  Range  Province,  is  characterized  by  high  temperatures  in 
the  lower  crust  and  upper  mantle.  This  region,  which  is  known 
to  extend  up  the  western  part  of  Canada,  is  also  a  poor  trans¬ 
mitter  of  long  period  surface  waves.  A  particularly  good 
example  of  this  phenomenum  can  be  seen  in  figure  11.  The 
Greenland  Sea  event  (table  1),  as  recorded  at  Grand  Saline, 
is  compared  to  the  same  event  recorded  at  Queen  Creek,  Arizona. 
This  station  (Fix  and  Sherwin,  1970)  is  located  in  the  heart 
of  the  Basin  and  Range  Province.  The  gains  of  the  seismographs 
at  the  two  stations  were  equalized  at  50  sec,  but  it  will  be 
seen  that  the  longer  period  surface  waves  are  virtually  miss¬ 
ing  at  the  Queen  Creek  site. 

To  the  east  the  continental  United  States  are  characterized 
by  the  Atlantic  and  Gulf  coastal  plains.  Furthermore,  the 
central  plains  are  known  to  have  relatively  low  temperatures 


in  the  lower  crust  and  upper  mantle  (Roy,  et  al.,  1968, 

Combs  and  Simmons,  1973)  which  provides  an  exceptionally  good 
path  for  Rayleigh  waves  with  periods  from  20-75  sec.  Our 
detection  station  in  north  Texas  is,  essentially,  at  the 
southern  end  ot  a  surface  wave  guide  which  passes  over  the 
north  pole  (Figure  12) . 

The  analog  plots  of  the  seismograms  used  in  the  study 
were  analysed,  after  they  had  ber  .  high-pass  filtered,  for 
the  presence  of  detectable  L  v  "ves.  Using  the  formula 


M  =  log  +  1.66  log  A  -  0*18,  which 

can  be  easily  derived  (Gutenberg,  1945),  values  of  Ms  were 
found  where  A  is  the  peak-to-peak  amplitude,  in  m/j,  of  the 
surface  wave  at  a  period  of  T=20  sec,  and  A  is  the  distance 
from  epicenter  to  detector,  in  degrees.  The  results  of  these 
computations  are  given  in  tables  1  and  3,  where  the  values  of 
Ms  are  compared  to  the  values  of  m^  determined  by  NCAA. 

Using  the  data  presented  in  this  study  an  interim  detec¬ 
tion  threshold  has  been  determined  for  our  Grand  Saline  record¬ 
ing  station.  For  events  travelling  within  the  polar  waveguide 
this  threshold  is  at  Ms  =  3*5,  for  a  single  vertical  seismometer, 
which  is  equivalent  to  a  threshold  of  m^  =  4*5  for  events  at 
depths  of  less  than  50  km. 
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FIGURE  CAPTIONS 


Response  of  the  seismograph  system. 

Rayleigh  wave  dispersion  information  recorded  in  north 
Texas  for  three  events  (see  table  1),  (a)  North  of 
Svalbard,  (b)  Greenland  Sea,  (c)  Norwegian  Sea. 

Rayleigh  wave  dispersion  information  recorded  in  north 
Texas  for  an  event  in  Sinkiang  Province  (see  table  1). 

Rayleigh  wave  dispersion  curves  for  the  four  events  in 
table  1  after  smoothing  (a)  North  of  Svalbard  (b)  Green¬ 
land  Sea  (c)  Norwegian  Sea  (d)  Sinkiang  Province. 

Composite  Rayleigh  wave  dispersion  curve  for  the  four 
events  in  table  1  in  conjunction  with  continental  disper¬ 
sion  curve  due  to  Oliver  (1262). 

Typical  dispersed  sine  wave,  or  chirp  filter,  derived 
from  the  dispersion  curve  shown  in  figure  4  (b) . 

Convolution  of  the  chirp  filter  derived  from  the  curve 
in  figure  5  with  each  of  the  four  events  given  in  table  1 
(a)  North  of  Svalbard,  (b)  Greenland  Sea  (c)  Norwegian 
Sea  (d)  Sinkiang  Province. 

Convolution  of  the  chirp  filter  derived  from  the  curve 
in  figure  5  with  three  of  the  events  given  in  table  3* 
events  from  very  different  locations. 


Convolution  of  the  chirp  filter  derived  from  the  curve  in 
figure  5  with  an  event  from  the  Hindu-Kush  (table  3)  in 
which  surface  waves  are  barely  discernible. 


Convolution  of  filter  and  surface  waves  showing  two,  or 
possibly,  three  surface  wave  arrivals  (a)  USSR  Mongolian 
border  on  17  May  1970,  =  4.5,  (b)  USSR  -  Mongolian 

border  on  23  May  1970,  xife  =  4.5.  (See  table  3  for  details). 
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11.  Comparison  of  the  surface  waves  from  the  Greenland 
Sea  event  (table  1)  recorded  at  (a)  Queen  Creek, 
Arizona,  (b)  Grand  Saline,  Texas. 

12.  A  polar  plot  of  the  world  centered  on  Grand  Saline. 
The  locations  of  the  epicenters  of  the  events  used 

in  this  study  are  shown;  the  alphanumeric  designation 
refers  to  tables  1  and  3. 
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TABLE  CAPTIONS 

Data  for  the  four  events  used  to  calculate  the 
composite  dispersion  curve  shown  in  figure  5. 
Coordinates  of  the  points  used  to  plot  the  curve 
in  figure  5. 

Data  for  the  other  events  used  in  this  study. 
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TABLE  2 


Group  Velocity 

Period 

3.9250 

75.5223 

3.8750 

69.2352 

3.8250 

60.8679 

3.7750 

51.8606 

3.7250 

45,6151 

3.6750 

41.9743 

3.6250 

39.5603 

3.5750 

37.4096 

3.5250 

35.5422 

3.4750 

34.0502 

3.4250 

32.0733 

3.3750 

30.5458 

3.3250 

29.3486 

3.2750 

27.7666 

3.2250 

26.1123 

3.1750 

24.6726 

3.1250 

23.4360 

3.0750 

21.3606 

3.0250 

19.2884 

2.9750 

18.6214 
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